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Abstract 

Background: Eliminating malaria requires vector control interventions that dramatically reduce adult mosquito 
population densities and survival rates. Indoor applications of insecticidal nets and sprays are effective against an 
important minority of mosquito species that rely heavily upon human blood and habitations for survival. However, 
complementary approaches are needed to tackle a broader diversity of less human-specialized vectors by killing 
them at other resource targets. 

Methods: Impacts of strategies that target insecticides to humans or animals can be rationalized in terms of 
biological coverage of blood resources, quantified as proportional coverage of all blood resources mosquito vectors 
utilize. Here, this concept is adapted to enable impact prediction for diverse vector control strategies based on 
measurements of utilization rates for any definable, targetable resource subset, even if that overall resource is not 
quantifiable. 

Results: The usefulness of this approach is illustrated by deriving utilization rate estimates for various blood, resting 
site, and sugar resource subsets from existing entomological survey data. Reported impacts of insecticidal nets 
upon human-feeding vectors, and insecticide-treated livestock upon animal-feeding vectors, are approximately 
consistent with model predictions based on measured utilization rates for those human and animal blood resource 
subsets. Utilization rates for artificial sugar baits compare well with blood resources, and are consistent with 
observed impact when insecticide is added. While existing data was used to indirectly measure utilization rates for 
a variety of resting site subsets, by comparison with measured rates of blood resource utilization in the same 
settings, current techniques for capturing resting mosquitoes underestimate this quantity, and reliance upon 
complex models with numerous input parameters may limit the applicability of this approach. 

Conclusions: While blood and sugar consumption can be readily quantified using existing methods for detecting 
natural markers or artificial tracers, improved techniques for labelling mosquitoes, or other arthropod pathogen 
vectors, will be required to assess vector control measures which target them when they utilize non-nutritional 
resources such as resting, oviposition, and mating sites. 
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Background 

While antiparasitic drugs and vaccines will be essential 
to the final stages of malaria elimination, their effective- 
ness as transmission control interventions will rely heav- 
ily upon unprecedented levels of vector control in highly 
endemic settings [1-4]. It will not be possible to elimin- 
ate malaria transmission from most of the tropics 
without developing scalable vector control intervention 
options which complement long-lasting insecticidal nets 
(LLINs) and indoor residual spraying (IRS) by targeting 
adult mosquitoes when they use resources other than 
human blood indoors, and indoor resting sites [5-9]. 
Mosquitoes usually also exploit non-human blood and hu- 
man blood sources outdoors, as well as sugar, outdoor 
resting sites, oviposition sites, and mating sites, so all of 
these other biological and environmental resources repre- 
sent alternative targets for vector control interventions [5] . 

Faced with such an array of resource target options, 
the challenge is to define exactly which of these inter- 
vention targets are optimal in each of the diverse vector- 
ial systems that exist [10], and to classify these settings 
into limited numbers of distinct categories where spe- 
cific intervention combinations maximize impact. Recent 
analyses indicate that the impact of vector control mea- 
sures targeting the blood hosts upon which mosquitoes 
depend can be rationalized in terms of measurements of 
the biological coverage of all available blood sources that 
is achieved, rather than merely high demographic cover- 
age of a targeted subset, such as humans while asleep 
indoors [9,11]. Blood resources are perhaps the best 
understood, and most readily quantified, of all the re- 
sources used by mosquitoes that could potentially be tar- 
geted with vector control measures [5]. However, many 
established or emerging vector control strategies only 
target a specific subset of the resource that can be read- 
ily identified and treated in the field. Examples include 
targeting insecticides to indoor resting sites only with 
IRS [12], and artificially introduced resting sites [13-17] 
or sugar baits [18-20]. Furthermore, the other available 
forms of these resources that cannot be targeted with in- 
secticides but compete with these subsets for the atten- 
tions of mosquitoes, such as naturally occurring sugar 
sources [21] or outdoor resting sites [22,23], are often 
impossible to identify or quantify [5]. It is therefore not 
possible to estimate biological coverage as a fraction of 
all available forms of that resource. Here, the concept of 
biological coverage is extended beyond blood resources, 
and adapted to enable impact prediction for more di- 
verse vector control strategies, based on direct measure- 
ments of coverage and utilization rates for definable, 
targetable subsets of less readily quantified resources 
that are equally important to mosquito survival, and 
therefore equally valid as potential targets for vector 
control interventions. 



Methods 

All symbols used are listed and defined in Table 1. 

Defining biological coverage based on the example of 
blood resources 

Biological coverage (C A>p ) of all the blood resources 
upon which mosquitoes rely, with long-lasting insecti- 
cidal nets (LLINs) or any other personal protection 
measure, has been previously defined as the proportion 
of all mosquito attacks upon all available hosts for which 
those hosts were covered with a protective (p) interven- 
tion at that time and place [11,24]: 

Ca, p =^ (1) 

where the total attack availabilities of the all hosts (A), 
and covered hosts at times and places when they are ac- 
tually protected (A p ), are defined kinetically [25,42] as 
the rates per night at which an individual host-seeking 
mosquito respectively encounters and attacks [26] either 
all hosts or all hosts that are protected at the time of the 
encounter and attack events. 

However, to allow simplified notation for generaliza- 
tion of this approach to a greater diversity of distinct re- 
sources, here the symbol A is replaced by R to denote 
the total kinetic availability rate of a specific given re- 
source, which may be specified as blood (v), resting sites 
(r), sugar (s) or any other resource (R e v, r, s ...). Also, 
the terms attack and protected which were previously 
used to define availability in kinetic terms for models of 
blood resource seeking and acquisition [26,27] are not 
entirely appropriate for non-blood resources, so these 
are replaced with more generally applicable terms at- 
tempt to utilize and covered, respectively. Furthermore, 
now that biological coverage has been defined to expli- 
citly consider only protection that is in place at the times 
and places when that resource is utilized, de facto cover- 
age (c) and protection (p) are equivalent to each other 
(c ~ p), so the former is used to simplify and harmonize 
the notation. Expressing Equation 1 in terms of this re- 
vised notation yields: 



where C v is the proportion of all mosquito attacks upon 
real (live vertebrate hosts) or perceived (artificial odor- 
baited traps, sometimes referred to as pseudo-hosts [27]) 
blood resources to which effective coverage with a vec- 
tor control intervention applies at that time and place, v 
is the total rate at which individual mosquitoes encoun- 
ter and attack all hosts and pseudo hosts, and v c is the 
total rate at which individual mosquitoes encounter and 
attack all hosts and pseudo hosts at times and places 
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Table 1 Parameter symbols and definitions 

Symbol Definition 

Previous formulations to predict impact based on of biological coverage of all blood resources 

A Availability of all blood hosts for attack, expressed as the rate at which they are collectively encountered and attacked per host-seeking 

mosquito per night [1 1,24-30] 

A hiP Availability of all protected (p) human (h) blood hosts for attack, expressed as the rate at which they are collectively encountered and 
attacked while protected by an LLIN or other prevention measure per host-seeking mosquito per night [1 1,24-30] 

C Ap Proportional coverage of all available blood resources that mosquito population utilizes (A) with a protective intervention (p) [1 1] 

Reformulation to predict impact based on coverage and utilization rates of resource subsets 

a R Utilization rate for an entire given resource (ft, which may be specified as blood (v), resting sites (r), sugar (5) or any other resource 

mosquitoes use), defined as the rate at which individual mosquitoes attempt to utilize all forms of that resource per gonotrophic cycle 

cxr x Utilization rate for a defined subset of a given resource (ft, which may be specified as blood (v), resting sites (r), sugar (s) or any other 
resource mosquitoes use) that can be identified and targeted with an intervention (x) in the field (ftj, defined as the rate at which 
individual mosquitoes attempt to utilize the subset per gonotrophic cycle 

<xr x Utilization rate for a defined subset (x) of a given resource (ft, which may be specified as blood M, resting sites (r), sugar (s) or any other 

resource mosquitoes use) that can be identified and targeted with an intervention (x) in the field during the subset of times (y) when it can 
be effectively protected by a given intervention (ft xy ), defined as the rate at which individual mosquitoes attempt to utilize that subset at 
times when it can be protected per gonotrophic cycle 

a/j x c Utilization rate for a defined subset of a given resource (ft, which may be specified as blood M, resting sites (r), sugar (5) or any other 

resource mosquitoes use) that has been identified, targeted (x) and covered (c) with an intervention during the subset of times (y) when it 
can be effectively protected by a given intervention (ft wc ), defined as the rate at which individual mosquitoes attempt to utilize that 
covered subset at times and places at which it can be protected per gonotrophic cycle 

<xr xj Utilization rates for a defined subset of a given resource (ft, which may be specified as blood (v), resting sites (r), sugar (5) or any other 

resource mosquitoes use) that has been identified, can be targeted with an intervention (x) and has been surveyed entomologically (z) in 
the field (R xz ), defined as the rate at which individual mosquitoes attempt to utilize that sample of that subset per gonotrophic cycle 

a,, Utilization rate for all available blood resources (v), defined as the rate at which individual mosquitoes attempt to utilize any source of 
blood per gonotrophic cycle 

a Vxz Utilization rate for a defined subset (x) of all blood resources (v) that has been identified and surveyed entomologically (z) in the field (v KZ ), 
defined as the rate at which individual mosquitoes attempt to utilize that sample of that blood source subset per gonotrophic cycle 

b Mean lifetime total number of bloodmeals acquired per emerging mosquito [26] 

6/ Mean mosquito biting rates experienced by individual livestock (/), defined as the number of bites per head per night 

B h Mean mosquito biting rates experienced by individual humans (/r), defined as the number of bites per person per night [26,30] 

C R Coverage of all available forms of a given resource (ft) with a vector control intervention 

Q, Coverage of all available forms of an identifiable, targetable subset (x) of a given resource (ft) with a vector control intervention 

C Vh Coverage of the human subset (x = h) of all available blood sources (v) while they are indoors (/) with long-lasting insecticidal nets (n) [1 1,26-30] 

E Emergence or recruitment rate of mosquitoes in a defined setting per night [24,26-29,31] 

h Humans 

/ Indoors 

Gonotrophic age, expressed as the number of gonotrophic cycles completed 

/ Livestock 

X t Relative availability of an individual mosquito traps (t) for attack by host-seeking mosquitoes attempting to utilize it as a source of blood, 
compared to a single unprotected human [27] 

M Absolute size of the mosquito population in a given setting, defined in terms of the number of individuals present 

m^ z Rate at which the mosquito population utilizes a defined, entomologically surveyed sample subset (z) of any identifiable and targetable 
subset (x) of a given resource (ft xz ), expressed as the number of utilization attempt events per night 

m„. Rate at which the mosquito population utilizes a defined, entomologically surveyed sample (z) of human (h) blood resources (v hz ), 
expressed as the number of utilization attempt events per night 

m„, Rate at which the mosquito population utilizes a defined, entomologically surveyed sample (z) of livestock (I) blood resources (V;J, 
expressed as the number of utilization attempt events per night. 

u Rc Mortality probability associated with exposure to an intervention-covered (c) form of a given resource (ft) through a single utilization 
attempt event 

A// Number of livestock (/) living in a defined setting [24,25,27-29] 

N h Number of humans {h) living in a defined setting [24,25,27-29] 
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Table 1 Parameter symbols and definitions (Continued) 

N hz Number of persons directly sampled by an entomological survey (z) of mosquitoes attacking human (h) hosts 

Nh,a Number of persons residing in all houses sampled by an entomological survey (O) of mosquitoes attacking human (h) hosts 

N, Number of mosquito traps (f) present in a defined setting [27] 

P n% Probability of a mosquito surviving all attempts to utilize intervention-covered forms of the targeted resource per gonotrophic cycle 
P Y Probability of a mosquito surviving all utilization attempt events for all resources per gonotrophic cycle [24,26,28] 

P Yr o Probability of a mosquito surviving all utilization attempt events for all resources per gonotrophic cycle in the absence of any intervention 
P f Probability of a mosquito surviving one full feeding cycle (/) [24,26-29,31] 

n Vhj Proportion of human (h) blood (v) host exposure to mosquito bites that occurs indoors (/) in the absence of any protective intervention 
[11,24,26,30,32-38]. 

Q Vc Proportion of all available bloodmeals (v) that originate from a specific livestock (/) host species subset [25] 

Q Vh Proportion of all available bloodmeals (v) that originate from the human (h) host species subset [1 1,22-27,30,39-41] 

ft The total availability of all forms of a given resource, which may be specified as blood ly), resting sites (r), sugar (5) or any other resource 

mosquitoes use, defined as the per night rate at which individual mosquitoes encounter and attempt to utilize that resource 

R x The total availability of a subset (x) of a given resource (ft which may be specified as blood (y), resting sites if), sugar (s) or any other 

resource mosquitoes use) that can be identified and targeted with an intervention, defined as the per night rate at which individua 
mosquitoes encounter and attempt to utilize that subset 

R xy The total availability of a subset (x) of given resource (ft which may be specified as blood (y), resting sites (r), sugar (5) or any other resource 
mosquitoes use) that can be identified and targeted with an intervention during the subset of times (y) when it can be effectively 
protected by that intervention, defined as the per night rate at which individual mosquitoes encounter and attempt to utilize that subset 
at times when it can be effectively covered with that intervention 

R Kyc The total availability of all intervention-covered (c) forms of a targetable subset (x) of given resource (ft which may be specified as blood (v), 
resting sites (r), sugar (5) or any other resource mosquitoes use) during the subset of times (y) when it can be effectively protected by that 
intervention, defined as the per night rate at which individual mosquitoes encounter and attempt to utilize the covered forms of that 
subset at times when it is effectively covered with that intervention. 

ft x-z The total availability of an entomologically surveyed sample (z) of a targetable subset (x) of given resource (ft which may be specified as 
blood M, resting sites W, sugar (s) or any other resource mosquitoes use), defined as the per night rate at which individual mosquitoes 
encounter and attempt to utilize it. 

r The total availability of all forms of resting sites, defined as the rate at which individual mosquitoes encounter and attempt to utilize resting 

sites per night 

s The total availability of all forms of sugar, defined as the rate at which individual mosquitoes encounter and attempt to utilize sugar per 

night 

t Mosquito traps 

T Mean number of nights individual mosquitoes spend resting and gestating indoors following a bloodmeal inside a house 

v The total availability of all forms of blood, defined as the rate at which individual mosquitoes encounter and attempt to utilize blood per 

night 

x A subset of a given resource that may be identified and targeted with a vector control intervention 

y A subset of a given resource that may be effectively covered with a vector control intervention at times and places when mosquitoes 

encounter and attempt to utilize it 

z A sample of a given resource that has been surveyed entomologically 

O Humans in a sampled set of households 



when they are effectively covered with a vector control 
intervention. 

In the case of interventions such as LLINs, which only 
protect humans while they use them indoors, biological 
coverage can be calculated as the product of the propor- 
tion of all bloodmeals (v) that originate from the human 
(h) host species subset (Q Vl/ ), the proportion of human 
exposure to mosquito bites that would otherwise occur 
indoors (z) without an LLIN {rc Vh A, and the proportional 
demographic coverage of humans, measured as the 



proportion of humans who reported using a net while 
asleep indoors the previous night (C Vji .) [11]: 

Cy = Qv^Vftj^,. (3) 

where all three terms on the right hand side of Equation 3 
are defined as sequentially nested fractions and sub- 
fractions of the total availability of all blood hosts (v) that 
are represented by humans (h), those humans while 
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indoors (/), and those humans while indoors and protected 
by coverage with LLIN use at the time (c): 



Qv h = Vh/V 



Cv hi = V h,i,c/Vh,i 



(4) 
(5) 
(6) 



Adapting the concept of biological coverage to rationalize 
vector control impact based on utilization rates of diverse 
resource targets 

Expressing Equation 2 in more general terms that may be 
applied to any given resource {R), rather than blood 
specifically (v), yields the following formula: 



C R 



Rc 
R 



(7) 



LLINs that directly kill mosquitoes when they encoun- 
ter and attack protected human blood sources are the 
best established [43] and easiest resource targets to 
conceptualize and model, but previous formulations pre- 
dicted their impact upon vector survival by assuming 
mosquitoes feed once and only once per gonotrophic 
cycle [26]. However, resting sites, oviposition sites, and 
even blood resources themselves, may be utilized more 
than once per gonotrophic cycle [44], while sugar 
sources may be used less than once [21,45]. To adapt 
the concept of biological coverage to more diverse re- 
source targets which are used more than once per gono- 
trophic cycle, the term resource utilization is defined as 
the mean rate at which mosquitoes utilize any given re- 
source (R) per gonotrophic cycle (a R ). This definition of 
resource utilization rate can be expressed mathematic- 
ally as the product of the duration of the gonotrophic 
cycle, expressed as nights per gonotrophic cycle (g), and 
the rate per night at which a mosquito population uti- 
lizes that resource (trig), divided by the size of the mos- 
quito population (M): 



a R 



g m R 
M 



(8) 



where g is expressed as units of nights, m R in units of 
utilization attempt events per night, and M as the number 
of individual adult mosquitoes present in the population. 
For any targetable, intervention-covered (c) proportion of 
that resource (R c ), the corresponding rate at which mos- 
quitoes encounter and attempt to utilize that covered 
fraction (ccr c ) , by definition, varies in proportion to the 
fraction of the kinetic availability of that resource that it 
represents: 



g m Rc g m R R c 



M 



M R 



(9) 



Hence the quotient of the availability or utilization 
rates for the total resource, divided into those for the 
intervention-covered fraction, are equivalent to bio- 
logical coverage of that resource: 



OCR c 
OCR 



m Rc 
m R 



Rc 
R 



(10) 



Most vector control strategies only target a specific 
subset (x) of the resource that they are delivered to, 
which is practically definable, identifiable, accessible, and 
treatable in the field. Similarly to resource coverage 
(Equation 10), the proportion of all available forms of a 
specific resource (R) that is accounted for by a given 
subset (R x ) of that resource (Qr x ) , can be defined and 
measured in terms of the rate at which mosquitoes en- 
counter and attempt to utilize it [26,27,42] by generaliz- 
ing Equation 4 for subsets of any possible resource, 
rather than specifying blood: 



R oc R 



OCR, 



(11) 



Similarly, the proportion of that subset (x) that is ef- 
fectively protected at times when it is utilized by mos- 
quitoes (y) can be expressed in terms of the proportion 
of resource utilization attempt events it accounts for in 
that resource subset: 



n Rx,y — a l*x,y I a K 



(12) 



Hence, the biological coverage of all forms of that resource 
(C R ) can be expressed more explicitiy than in Equation 10 as 
the product of the proportion of that resource represented 
by that subset (Qr ) , the proportion of utilization attempt 
events for that subset to which protection effectively and 
conditionally applies (hr xt ) and measured intervention 
coverage of that resource subset at times and places when it 
may be effectively protected (Cr xj ): 



ocr x 



Cr — Q Rx JiR xy CR xy — — — Cr, 



ocr x 



ocr 



OCR 



(13) 



Note that CLr is the utilization rate of the covered 
fraction of the targeted resource subset, equivalent to 
a Rc in Equation 10, because all covered forms of the re- 
source occur amongst the intervention-targeted subset 
of that resource (x) at the times and places at which they 
mosquitoes actually attempt to utilize them iy): 



otR xyc = cxr c because R C ^R. 



(14) 



Predicting intervention impact based on resource subset 
coverage and utilization rates 

Interventions targeting adult mosquitoes may have quite 
complex modes of action, repelling mosquitoes away 
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from humans [46] or contaminating them with agents that 
affect their longevity [47,48], competence [47,48] or fe- 
cundity [49,50]. Biological agents may be transmitted hori- 
zontally or vertically through the population [47,48,50], 
while coverage amplification of chemicals may be achieved 
by mosquito-mediated transfer between resources [49]. 
Regardless of their complexity, these diverse strategies can 
all be enhanced by maximizing biological coverage of the 
resource targeted to ensure maximum contact with the 
mosquito population, and this is a critically important de- 
terminant of success in its own right. Previous formula- 
tions describing mosquito survival and mortality as a 
function of exposure to LLINs or IRS [26] are therefore 
adapted and simplified as follows to allow for a range of 
utilization rates ranging from zero to several times per 
gonotrophic cycle, rather than the previously assumed 
utilization rate of once per gonotrophic cycle for all blood 
resources (a v =l). All predictions of impact upon mosquito 
survival, and the entomologic inoculation rates they medi- 
ate, were implemented and parameterized exactly as pre- 
viously described [51], except that Equation 14 of the 
original formulation [26] was adapted to incorporate the 
mortality risks of utilizing all covered and uncovered re- 
sources in a more generally applicable manner: 

Py = PyfiPa llc (15) 

where P y is the probability of surviving all utilization at- 
tempt events for all resources per gonotrophic cycle, P YiU 
is the probability of surviving all utilization attempt 
events for all resources per gonotrophic cycle in the ab- 
sence of any intervention, and P ai!c is the probability of 
surviving all attempts to utilize intervention-covered 
forms of the targeted resources per gonotrophic cycle. 
The probability of surviving all attempts to utilize 
intervention-covered forms of the targeted resource per 
gonotrophic cycle [P aRc ) can be calculated as an expo- 
nential decay function of the product of the mortality 
probability associated with exposure to a covered form 
of the resource through a single utilization attempt event 
(fi R ) , and the mean utilization rate for all covered forms 
of that resource (or c ): 

Pa Bc = e-W«° (16) 

By substituting rearranged forms of Equation 10 and 
then Equation 13 into Equation 16, a solution with two 
field measurable parameters for the targetable, quantifi- 
able, surveyable, subset is derived: 

P ag = e -PXc a * C >> = e ~l i Rc a «x.y C B x ,y ^7) 

It is therefore not necessary to know the proportion of 
that total resource which the targeted subset represents, 
or the coverage (Cr) or utilization rate (or) for all avail- 
able forms of a resource. Impact can be predicted 



directly so long as the coverage of the targeted subset it- 
self (Cj^J, and utilization rates for that subset under con- 
ditions that enable the intervention to protect it against 
safe utilization by the mosquito (oCr ), can be measured. 

This approach to predicting the survival probability as- 
sumes that utilization attempt events are randomly, and 
independently, distributed across all resource units and 
mosquitoes. Specifically, the number of times one mos- 
quito utilizes a resource (or resource subset) in one gono- 
trophic cycle is assumed to be a non-negative integer 
valued random variable (0, 1, 2, 3...) since the mosquito 
may not necessarily use the resource or, alternatively, may 
access it multiple times. Hence, the utilization rate of 
these resources should be understood as an expected 
value depending on random events that may be expressed 
as a mean. This is clearly not the case in relation to obli- 
gate utilization of blood from one of all available blood re- 
sources (R = v). Each mosquito must utilize one of these 
available resources to complete the gonotrophic cycle, so 
complete coverage of all blood resources (Cr = 1) that are 
utilized at a mean rate of once per gonotrophic cycle 
((Xr - 1) with an insecticide which induces comprehensive 
fatality (u R = l) would deterministically result in reduction 
of survival probability to zero (P asc = <*rCr = 0), ra- 
ther than merely reduced to the minor proportion of mos- 
quitoes that are inaccurately assumed by Equation 16 to 
have completed a gonotrophic cycle without taking any 
bloodmeal. However, for a covered subset of a resource 
(Equation 17), rather than for all available forms of that 
resource (Equation 16), it is realistic to assume that the 
number of utilization attempt events per gonotrophic 
cycle is a random variable for individual mosquitoes and 
utilization rates per gonotrophic cycle are expected values 
(expressed as a mean), even for obligate blood resource 
utilization behaviours. 

Measuring utilization rates for subsets of undefined 
resources by comparison with those for quantifiable 
blood resources 

Adult mosquitoes use many distinct resources during 
their lifetimes, including several that they need afresh 
every time they complete a gonotrophic cycle: blood, 
resting sites, and oviposition sites. Most of these re- 
sources are difficult to quantify directly, so the same is 
true of the rates at which mosquitoes utilize them, 
thereby making contact with them. However, measure- 
ments of feeding rates upon humans or livestock allow 
ready quantification of absolute mosquito population 
size or recruitment rate [31]. This is because the size of 
these mammalian host populations can be conveniently 
measured by direct census, and blood acquisition occurs 
at a measurable rate per host [52] with a measurable 
probability for a given blood host species [22,23]. Also, 
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blood acquisition usually occurs at a utilization rate of 
only once per gonotrophic cycle (a v ~ 1) where sugar avail- 
ability is not limiting [45], except for the first gonotrophic 
cycle where two bloodmeals may be required [53,54]. For 
example, the emergence or recruitment rate of mosqui- 
toes (£) in a given setting can be calculated as function of 
the measured mean biting rate experienced by individual 
humans (£>/,), the number of humans living there (TV/,), the 
proportion of bloodmeals obtained from humans (Q n ) , 
and survival probability per feeding cycle (Pj) [31]: 



B h N h 



(18) 



where the mean lifetime number of bloodmeals per 
emerging mosquito (b) is calculated as the sum of the 
probabilities of surviving to all plausible gonotrophic 
ages, expressed as the number of gonotrophic cycles 
completed (/') [26,31]: 



/ 



(19) 



Similarly, for a very zoophagic (predominantly animal- 
feeding) vector with a strong preference for a known, ac- 
cessible, manageable non-human host species such as 
cattle, goats, sheep, pigs or other livestock (/), it may be 
easier to accurately measure biting rates upon such live- 
stock (5;) so the equivalent calculation can be made if 
the proportion of blood obtained from that host species 
(Q V( ) , and the population size of that host species (A//) 
can be determined: 



BiN, 
Q Vl b 



(20) 



The key to applying Equations 18 and 20 to estimate ab- 
solute vector population sizes is the assumption that the 
fraction of all available sources of blood that each entomo- 
logically surveyed host represents [Q v /Nf, or Q v /AT/) 
can be readily estimated by host census and bloodmeal 
identification from a sample of resting, blood-fed mosqui- 
toes, so it is not necessary to directly detect all biting 
events on all hosts. Generalizing this principle, the emer- 
gence rate of mosquitoes (£) can be estimated based on 
the rate at which mosquitoes are trapped or observed util- 
izing (»Jr, z ) a surveyed sample subset (z) of any targetable 
subset (x) of a given resource {R x>z ) if the proportion of all 
available forms of that resource which that surveyed sub- 
set represents (R XrZ /R), and the rate at which individual 
mosquitoes utilize all available forms of that resource 
per gonotrophic cycle (an), are both known. Note also 
that equation 18 and 20 both also implicitly include a 
term in the denominator for the utilization rate of all 
blood sources, that was negated by assumed a value ap- 
proximating unity (a v ~ 1) but can be explicitly 



reintroduced for the purposes of generalization. Substi- 
tuting m Rxz for B h or B b R XiZ /R for Q v jN h or Q n /Ni, 
and a R for a v in Equations 18 and 20, respectively, yields 
the following general formula: 



m R x , 



R 

Ry > 



/ (a R b) 



(21) 



Blood resources can be readily identified as discrete 
units and their total numbers can be quantified by head- 
count census. However, units of sugar, resting site, ovi- 
position site, and mating site resources are difficult to 
define unambiguously, except where these are intro- 
duced artificially (sugar baits, houses, boxes, pots, bar- 
rier screens, water containers, or swarming markers), 
and the total quantity of these resources available in the 
environment is even more difficult, if not impossible, to 
ascertain. Therefore, it is not obvious how R XlZ /R can be 
estimated for these resources with existing field survey 
methods. However, this is not necessary to know, because 
intervention impact can be conveniently rationalized in 
terms of coverage and utilization rates for definable, tar- 
getable subsets of those resources (Equation 17), and it is 
possible to calculate their relative rates of utilization com- 
pared with those for blood from a known proportion of all 
available blood resources. Here the quantifiable total 
blood resource, and a surveyed sample (z) of an identifi- 
able subset (x) of that blood resource, is specified (R = v) 
and distinguished from equivalent terms for other re- 
sources, such as resting sites (R = r) or sugar (R = s), with 
the specific terms v and v X:Z . Also, the per gonotrophic 
rates at which mosquitoes attempt to utilize (a) all blood 
resources (v) or a distinct, identifiable subset of blood re- 
sources (v x ), as well as the per night rate at which 
utilization attempt events occur (w) at a surveyed sample 
(z) of that blood resource subset (v x ), are distinguished 
from those for other resources with the specific terms a v , 
a Vx , and m v , respectively. Given that resource utilization 
rate per gonotrophic cycle is proportional to the total rate 
at which utilization attempt events occur in the overall 
population (Equations 8, 9, 10 and 11), the relative rate of 
utilization of such a resource subset compared to all blood 
resources can be expressed by dividing Equation 21, which 
is specified for a given non-blood resource subset and a 
surveyed sample thereof (R x /R xz ), by an equivalent formu- 
lation specified for all blood resources, and a surveyed 
sample of hosts from a subset for which bloodmeals re- 
covered from the midguts of recently fed specimens can 
be identified and distinguished from other sources (v/v xz ), 
and rearranging: 



<xr x _ m R> , _ m Rxz R x v x 
oc v 



(22) 
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Note that the emergence and mean longevity terms 
cancel each other out so that estimates of these parame- 
ters are not required to estimate the relative rate of 
utilization of a resource compared with blood, as de- 
scribed below. 

The most obvious vertebrate blood resource subsets 
{v x ) which are readily surveyed, including detection of 
blood feeding events and identification of blood source 
in specimens of fed mosquitoes, are humans (x = h) and 
livestock (x = I) [52]. The proportion of all available 
blood resources sampled by the field survey can be 
quantified as the product of the proportion of blood- 
meals obtained from humans (Q V J or livestock (Q Vl ) 
and the number of humans (N hiZ ) or cattle (N lrZ ) sampled 
by the host attack survey, divided by the total number of 
humans (N h ) or livestock (TV/) present: 



V 



N h ,zQ Vl 



or 



v 



Ni, z Q v , 



(23) 



Fortunately, while gonotrophic discordance beyond 
the first feeding cycle does occur in Anopheles, it is un- 
usual and can be quantified [45]. In most cases, it is 
therefore reasonable to explicitly assume that mosqui- 
toes predominantly utilize blood approximately once per 
feeding cycle (a v ~ 1), so substituting the host-specified 
( v x = v h or V/) formula of Equation 23 into Equation 22, 
and replacing a v with unity, yields a solution for ccr x for 
which all the terms are measurable in the field: 



OCR, 



m Rxz R x N Kz Q v 
m n R x . z N h 



m Rx R x N Lz Q Vi 
m Vl R x , z N t 



(24a) 



(24b) 



where R X JR X is the proportion of all available forms of 
the targeted non-blood resource subset that was sur- 
veyed entomologically to measure the rate per night at 
which the entire mosquito population attempts to utilize 
it (mR xz ), where N hiZ /N h and N l>z INi are the proportions 
of all humans or livestock that were respectively surveyed 
to measure the rate at which mosquitoes attempted to 
utilize their blood, and where Q Vh and Q v are the propor- 
tions of bloodmeals the vector population obtains from all 
available humans and livestock, respectively. 

Where two resources co-occur and overlap completely 
with each other, specifically the example of resting sites 
(R = r,) and human blood indoors within houses (vj,,,-)> 
the proportion of each resource subset that is sampled is 
no longer required because these cancel each other out. 
All that is required is an estimate of the number of per- 
sons or person nights sampled by the host attack survey 



{N hiZ ), and the total number people staying in those sam- 
pled houses {N^n), or even their ratio, which is commonly 
referred to as the mean number of occupants per house 



ocr x 



m Rxz N h . z Q Vi 
m v ,N h n 



(25a) 



In some experiments, however, resting events follow- 
ing more than one bloodmeal are represented in surveys 
of resting sites because those events may last two or 
more days. Recent standardized trials to compare vari- 
ous techniques for catching host-seeking and resting 
mosquitoes [55-57] placed these in or immediately out- 
side of different houses within a defined sampling frame 
each night, so that the former would not compete with 
the latter by trapping mosquitoes before they can feed 
and rest. Also, spray catches must be spaced by intervals 
of several days to allow residual pyrethrum to dissipate. 
In both cases, mosquitoes gestating over two or more 
preceding nights are allowed to accumulate from multiple 
nights of blood feeding (r) in the surveyed houses, and 
this must be accounted for when estimating utilization 
rates: 



CtR, 



m Rxz N h . z Q Vfi 



(25b) 



Literature review and utilization rate estimate extraction 

Studies, or sets of studies, were identified which 
presented sufficient parameter estimate data for 
utilization rates of specific, intervention-targetable re- 
source subsets to be calculated for specific malaria 
vector species in specific, distinct locations. In 
addition to the authors' archives of literature and un- 
published data, the Pubmed database was also queried 
with the search term 'Anopheles AND ((pyrethrum 
spray OR aspirator) OR (insecticide AND (cattle OR 
livestock)) OR odour-baited OR sugar)'. For utilization 
of livestock blood and sugar, consideration was lim- 
ited to studies in settings where trials of insecticide- 
treated livestock or sugar baits, respectively, have 
been either implemented or specifically suggested. To 
avoid cluttering of the second figure presented in the 
Results section, consideration of studies enabling esti- 
mation of indoor resting site utilization was restricted 
to recent unpublished studies of our own and those 
published in the last decade. Where results for a spe- 
cies complex or group were reported, these are attrib- 
uted to the most common sibling species identified in 
that population. 

Utilization rates for blood from humans while indoors, 
when they can be protected with LLINs, was calculated 
as the product of the proportion of human exposure to 
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mosquito bites occurring indoors (n Vhl ) and the human 
blood index (Q v .) , by assuming a single bloodmeal per 
gonotrophic cycle: 



n vuQv, 



(26) 



Where local estimates for the proportion of blood- 
meals obtained from humans (Q v;j ) were not available 
for the vector in question, the median values for that 
species from a previous review [39] were applied. Where 
direct estimates of the proportion of bloodmeals ob- 
tained from cattle and other treated livestock (Q VJ ) were 
not available, utilization of blood from other non-human 
sources was assumed to be negligible, so that this quan- 
tity could be calculated as the complement of the pro- 
portion obtained from humans (Q„J: 



(27) 



Utilization rates for odour-baited traps are calculated 
by assuming that the probability of a mosquito attacking 
a trap, rather than a natural host, per gonotrophic cycle 
is equivalent to the proportion of all available human 
hosts, animal hosts, and pseudo-hosts, that they repre- 
sent [27]. This is calculated on the basis of the ratio of 
traps to people (N t /N h ), the relative availability of those 
traps (X t ), and the proportion of bloodmeals obtained 
from non-human hosts (l-Q v J [27]: 



(l+((l-Q Vh )/Q Vh ) + (X t N t /N h )) 



(28) 



Utilization rates for sugar resource subsets (a Sx ) were 
calculated as follows, based on direct field measure- 
ments of utilization rates per day for dye-labelled sugar 
baits (a* ) [18-20], and an assumed mean gonotrophic 
cycle duration of three days (g = 3): 



a sg 



(29) 



Utilization rates for resting site subsets (a r J, such as 
the insides of houses or artificial shelters and netting 
barriers placed in or around them, were estimated using 
Equation 25a or 25b based on the quotient of the mean 
rate at which mosquitoes were caught resting inside a 
sample of them {m rXiZ ) adjusted, where necessary, for an 
assumed indoor resting period of 2 days (r = 2), divided 
by the rate at which they were caught attacking human 
hosts indoors (m Vhl ), the mean number of exposed occu- 
pants per house or room {N hin IN hiZ ) adjusted for reported 
usage rates of LLINs (C VKtn ) which were assumed to con- 
fer complete protection, the proportion of bloodmeals ob- 
tained from humans (Q V(I ) , and the proportion of human 
bloodmeals obtained indoors {n Vhi ) , with the latter two 



parameters usually assumed from mean literature values 
reported for that species [32,39]: 



m rxz n Vhi N h:Z Q Vh 



(30) 



Results 

Dependence of impact upon utilization rates of 
intervention-targeted resource subsets 

Figure 1 illustrates how the impact of a vector control 
intervention with a mosquito-toxic active ingredient de- 
pends on the utilization rate of the resource subset to 
which it is targeted, and on the pre-existing level of trans- 
mission mediated by the vector. Utilization rates exceed- 
ing, or at least approaching, one event per gonotrophic 
cycle are required to achieve useful reductions of intense, 
saturating transmission [4,58] mediated by human-feeding 
mosquitoes such as Anopheles arabiensis and Anopheles 
gambiae. Less efficient vectors that predominantly feed 
upon animals, such as Anopheles cuiicifacies, cause far 
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Figure 1 The predicted relationship between the impact of 
vector control upon malaria transmission and the rate at which 
mosquitoes utilize the resource subset targeted by the 
intervention Or x . Anopheles epiroticus, An. arabiensis, and An. gambiae 
were chosen as examples of vector species that feed primarily upon 
animals, on both humans and animals, or primarily upon humans. Zero 
coverage of the human blood resource with LLINs (C V)i = 0) and high 
biological coverage of the resource subset which is targeted with an 
intervention (Q, = 0.8) that induces high mortality among 
mosquitoes that utilize covered forms of that resource {u R = 0.8) 
were assumed. Apart from the adaptation of the original formulation 
to capture intervention impact upon vector survival as a generally 
applicable function of target resource subset coverage and utilization 
rate (Equation 1 5 to 1 7), all predictions were made by parameterizing 
and executing deterministic models of malaria transmission and 
control, exactly as previously described [26,51,62]. 
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lower levels of transmission and morbidity burden [39] 
that can respond more sensitively to effective transmission 
control because baseline exposure levels are not sufficient 
to saturate the human population with parasite infections 
[24,58-61]. Very valuable impact upon transmission by 
such zoophagic mosquitoes may even be achieved by tar- 
geting resource subsets that are used as rarely as once 
every three feeding cycles because most malaria is trans- 
mitted by mosquitoes that are at least four gonotrophic 
cycles old [31,54]. However, unless a resource subset is 
used at least once in every five gonotrophic cycles, it is un- 
likely to be a useful target for suppressing mosquito sur- 
vival and population density through vector control, even 
for weak zoophagic vectors that rarely feed upon humans 
and mediate modest, manageable levels of transmission. 

Field estimates of utilization rates for defined, targetable 
resource subsets 

Figure 2 illustrates the range of utilization rate estimates 
that could be extracted from the literature for a diversity 
of resource subsets used by mosquitoes. The predictions 
of Figure 1, and the generally high utilization rates of 



human blood while people are indoors (Figure 2), are 
approximately consistent with the documented impacts 
of LLINs upon malaria transmission [43]. It is also not- 
able that the lowest estimates for utilization of human 
blood indoors (Figure 2) are for An. ambiensis, Anoph- 
eles farauti, Anopheles darlingi, and Anopheles nunezto- 
vari, all of which have been known to persist and 
dominate residual transmission systems following scale- 
up of LLINs or IRS [6-8,33]. Insecticide-treated livestock 
also appear to be as promising a target for zoophagic 
vectors as LLINs are for anthropophagic (predominantly 
human-feeding) vectors, and the utilization rates esti- 
mated for Anopheles culifacies and Anopheles stephensi 
are consistent (Figure 2) with the proven impact of this 
approach upon malaria transmission by these species 
[63]. Odour-baited traps act as pseudo-hosts by mimick- 
ing, and even surpassing [64], the taxis stimuli of normal 
blood sources for mosquitoes, luring them to fatal trap 
devices or insecticide-treated surfaces [27]. While these 
can achieve useful utilization rates if their attractiveness 
and positioning can be optimized, these estimates are 
somewhat lower than for insecticidal nets and livestock 
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Figure 2 Estimated utilization rates for a variety of blood, resting site and sugar resource subsets (a R J based on existing field data. 

See Additional file 1 for details of literature, parameter values and equations used to calculate the estimates presented. 
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treatments directed at the best-matched vector species 
(Figure 2) because these devices must compete with nat- 
ural hosts that, therefore, constitute an inevitable gap in 
biological coverage [27]. 

Apart from blood, the other important nutrition 
source that facilitates mosquito survival and malaria 
transmission is plant-derived sugar [21,65,66]. Estimated 
minimum utilization rates for dye-labelled sugar baits by 
Anopheles claviger [19], Anopheles sergenti [18], and An. 
gambiae [20], in three distinct settings appear to be at 
least comparable with those for utilization of human 
blood indoors by very anthropophagic vectors, and for 
utilization of animal blood by zoophagic vectors 
(Figure 2). The impressive impacts upon all three of 
these vector populations that have been achieved by 
adding insecticides to such sugar baits [18-20] are there- 
fore consistent with the predictions outlined in Figure 1, 
as well as other recent modelling analyses [67]. Given 
the widespread dependence of mosquitoes upon sugar 
[21,65,66], many important vector populations probably 
use it at similarly high rates, especially when infected 
with malaria parasites [68]. Mosquitoes should therefore 
be at least as amenable to control with this approach as 
anthropophagic vectors are to LLINs [43], and as zoo- 
phagic vectors are to insecticide-treated livestock [63]. 

Utilization rate estimates for indoor resting sites 
(Figure 2) are generally lower than that required to ex- 
plain (Figure 1) the often massive impact of IRS [13] and 
insecticidal wall linings [69] on many target vector spe- 
cies. However, this is not entirely surprising because 
even the best techniques for capturing mosquitoes rest- 
ing indoors, such as pyrethrum spray catch [52] and 
backpack aspirators [70], are known to consistently 
under-sample them. Nevertheless, the estimates of rest- 
ing site utilization presented in Figure 2 are clearly use- 
ful for comparison of different potential resting site 
targets, confirming that comprehensive spraying or lin- 
ing of entire rooms and houses is probably superior to 
targeted treatment of pots, boxes (some of which were 
baited with host odours) or screening barriers, all of 
which were placed indoors for endophilic vectors or out- 
doors for exophilic ones [13-17]. However, perhaps the 
most important observation in relation to these esti- 
mates of resting site subset utilization rates is that they 
rely on upon quite complicated calculations requiring at 
least five distinct input parameters (Equations 25a, 25b 
and 30), many of which have to be assumed based on 
best guesses or literature values from a different setting 
(Additional file 1). In fact, none of the estimates pre- 
sented in Figure 2 are based entirely upon local esti- 
mates for all the input parameters (Additional file 1), 
and are therefore not entirely independent of each other 
or representative of the full range of values for any of 
the vector species described. 



Discussion 

Defining and surveying targetable resource subsets 

Estimating coverage and utilization of a resource subset 
primarily depends upon defining it in a quantifiable 
manner that can be readily surveyed and targeted, or ar- 
tificially created in the field. The most obvious and fa- 
miliar examples are the human populations targeted for 
universal coverage with LLINs to protect the blood re- 
source they represent to mosquitoes [28,43,71]. While 
wild animals are difficult to survey or deliver interven- 
tions to, livestock represent blood resources that can be 
readily quantified and targeted with interventions [63]. It 
is even easier to track numbers and functionality of arti- 
ficial odour-baited traps, which mimic and compete with 
natural blood sources for the attentions of host-seeking 
mosquitoes, so their potential impact can also be pre- 
dicted as a function of biological coverage of all available 
hosts and pseudo-hosts [27]. 

The subset of all resting sites represented by the inner 
surfaces of human dwellings (walls, ceilings and even 
furniture) are the defined target for IRS [12,71], as well 
as insecticide-treated tents [72], shelters [73,74], and 
wall linings [69], so coverage can be quantified as the 
proportion of residential structures treated. On the other 
hand, artificially introduced pots, boxes, curtains, linings 
or screening barriers compete with natural resting sites 
[13-17]. Mosquitoes can be captured relatively efficiently 
on these well-defined, convenient, standardized surfaces, 
so it has been suggested that these could also be treated 
with toxic insecticides to improve control efficiency 
[13-17]. While it remains difficult to consistently identify 
and define sugar, oviposition site or mating site re- 
sources [5], recent progress with observational [75,76], 
trapping [77], tracing and labelling [18-20,22,23,78-81] 
methods for mosquitoes is encouraging. 

Adapting entomological survey techniques to measure 
resource utilization rates 

Comparing the range of utilization rates described in 
Figure 2 with the predictions of potential impact illus- 
trated in Figure 1 confirms that, despite their known 
limitations [52,82], existing entomological field methods 
may be very useful for designing and evaluating a wide 
diversity of vector control products [83]. Both blood and 
sugar meals can be readily identified using a variety of 
naturally-occurring markers and artificially added tracers 
[18-20,22,23,78-81], thus enabling very direct, robust 
measurement of label uptake as a function of time or 
age. Utilization rates can therefore be estimated directly 
for subsets of these naturally occurring resources 
(Equations 26, 27, and 29) or indirectly for artificially intro- 
duced subsets such as odour-baited traps (Equation 28). 

Utilization rates for resting site subsets (Equation 30), 
or indeed any other non-blood resource (Equation 25a 
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and 25b), can also be estimated indirectly by calibrating 
against measurable utilization rates for quantifiable, pre- 
ferred blood sources. However, the complexity of these 
models, and their reliance upon local measurements of 
several entomological input parameters, all of which 
have limited precision and accuracy, may well limit 
broader application of this approach beyond the crude 
application to existing data presented in Figure 2. Recent 
attempts to rejuvenate and improve existing entomo- 
logical survey methodology for detecting resource uti- 
lization attempt events with electrified grids [84-86], 
sticky traps [77,87], mechanized aspirators [70], and high 
resolution cameras [75,76], should enable improved sen- 
sitivity of utilization attempt event detection at surveyed 
samples of resource subset targets. However, while such 
technical advances may well address the inaccuracies of 
attempts to estimate utilization rates for subsets of rest- 
ing sites or other non-nutritional resources by improving 
event detection sensitivity, they are unlikely to improve 
their precision because considerable uncertainty arises 
from the need for relatively complex models (Equations 
25a, 25b and 30) that require correspondingly numerous 
measurements of input parameters. 

Fortunately, a wide range of more sensitive chemical, 
biochemical, genetic and biological markers, that could be 
applied to labelling mosquitoes when they use these other 
resources, are now available [82] but these remain to be 
fully exploited. In fact, field studies using artificial tracers 
to label of both mosquitoes feeding upon sugar [18-20] 
and sand flies feeding upon rodent blood [88,89], in which 
addition of insecticide removed almost all labelled insects 
from these vector populations, clearly demonstrate the 
validity of this strategy as a means to estimate biological 
coverage or utilization rates. The major advantage of la- 
belling mosquitoes when they utilize a resource subset, 
rather than trapping or observing them, is that the mea- 
sured proportions of marked insects can be readily ana- 
lyzed with robust off-the-shelf statistical methods for 
binary outcomes, and are relatively precise because they 
have a nominator and denominator which both vary in 
proportion to population size or event detection sensitiv- 
ity. The largest caveat to this approach is that essentially 
all targeted forms of that resource must be labelled on 
geographic scales large enough to negate the effect that 
mosquito dispersal has upon measurements of label up- 
take: immigration of unlabelled mosquitoes into the study 
area will increase the denominator while emigration will 
reduce the nominator, so that true local utilization rates 
will be systematically underestimated [52,90,91]. However, 
this phenomenon could also be exploited to great advan- 
tage if multiple distinct labels for various treatment arms 
were used to measure, and correct for, the effects of mos- 
quito upon impact distribution in large-scale trials of vec- 
tor control interventions [92] . 



The conceptual framework and entomological meas- 
urement priorities outlined here should be readily and 
directly applicable to almost any population of mosqui- 
toes, vectors or other pest. It should therefore be pos- 
sible to simultaneously tackle multiple vectors with 
integrated vector management programmes [93] that 
prioritize interventions based on simultaneous, com- 
parative field assessment of respective utilization rates 
for each potential target species. Recent demonstrations 
of the usefullness of dummy bait products containing 
appropriate labels but no insecticide [18-20,88,89] illus- 
trate how cost-effective, robust measurements of utiliza- 
tion rates could be used to select and optimize available 
technologies for immediate use or new prototypes for 
development. 

Conclusions 

The concept of biological coverage can be extended to en- 
able prediction of intervention impact for diverse vector 
control strategies based on estimated utilization rates for 
any definable, targetable resource subset. Indeed the ap- 
plicability of this approach has been demonstrated here 
using existing entomological measurement methods to ra- 
tionalize the observed impacts of LLINs, insecticide- 
treated livestock, and attractive toxic sugar baits upon 
malaria vectors. The development of improved and diver- 
sified technologies for controlling transmission of malaria, 
as well as a diversity of other vector-borne pathogens, 
could therefore be accelerated, rationalized and stream- 
lined based on field measurements of the rates at which 
mosquitoes utilize targetable biological resource subsets. 

While blood and sugar consumption can be readily 
quantified using existing methods for detecting natural 
markers or artificial tracers, improved techniques for la- 
belling mosquitoes will be required to assess and optimize 
vector control measures which target them when they 
utilize resting, oviposition and mating sites. All mosquito 
species need sugar, resting sites, oviposition sites, and 
mating sites, as indeed do most arthropods of medical and 
veterinary importance. These resources are therefore im- 
portant potential targets for the new or improved vector 
control methods that are clearly needed to eliminate mal- 
aria, and also a variety of other vector-borne pathogens. 
To enable comparative assessment of all potential re- 
source subset targets, including sites which mosquitoes 
rest, oviposit or mate at, existing tracer technologies need 
be adapted to enable reliable, non-toxic, non-disruptive la- 
belling of mosquitoes when they utilize these non- 
nutritional resource subset targets. 

Additional file 



Additional file 1: Data sources, publications and calculations 
supporting Figure 2. 
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